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The spread of disease through a physical-contact network and the spread of information about the disease on 
a communication network are two intimately related dynamical processes. We investigate the asymmetrical 
interplay between the two types of spreading dynamics, each occurring on its own layer, by focusing on 
the two fundamental quantities underlying any spreading process: epidemic threshold and the final 
infection ratio. We find that an epidemic outbreak on the contact layer can induce an outbreak on the 
communication layer, and information spreading can effectively raise the epidemic threshold. When 
structural correlation exists between the two layers, the information threshold remains unchanged but the 
epidemic threshold can be enhanced, making the contact layer more resilient to epidemic outbreak. We 
develop a physical theory to understand the intricate interplay between the two types of spreading dynamics. 

Epidemic spreading 16 and information diffusion 710 are two fundamental types of dynamical processes on 
complex networks. While traditionally these processes have been studied independently, in real-world 
situations there is always coupling or interaction between them. For example, whether large-scale outbreak 
of a disease can actually occur depends on the spread of information about the disease. In particular, when the 
disease begins to spread initially, individuals can become aware of the occurrence of the disease in their neighbor- 
hoods and consequently take preventive measures to protect themselves. As a result, the extent of the disease 
spreading can be significantly reduced 1113 . A recent example is the wide spread of severe acute respiratory 
syndrome (SARS) in China in 2003, where many people took simple but effective preventive measures (e.g., 
by wearing face masks or staying at home) after becoming aware of the disease, even before it has reached their 
neighborhoods 14 . To understand how information spreading can mitigate epidemic outbreaks, and more broadly, 
the interplay between the two types of spreading dynamics has led to a new direction of research in complex 
network science 15 . 

A pioneering step in this direction was taken by Funk et at, who presented an epidemiological model that takes 
into account the spread of awareness about the disease 16,17 . Due to information diffusion, in a well-mixed 
population, the size of the epidemic outbreak can be reduced markedly. However, the epidemic threshold can 
be enhanced only when the awareness is sufficiently strong so as to modify the key parameters associated with the 
spreading dynamics such as the infection and recovery rates. A reasonable setting to investigate the complicated 
interplay between epidemic spreading and information diffusion is to assume two interacting network layers of of 
identical set of nodes, one for each type of spreading dynamics. Due to the difference in the epidemic and 
information spreading processes, the connection patterns in the two layers can in general be quite distinct. 
For the special case where the two-layer overlay networks are highly correlated in the sense that they have 
completely overlapping links and high clustering coefficient, a locally spreading awareness triggered by the 
disease spreading can raise the threshold even when the parameters in the epidemic spreading dynamics remain 
unchanged 16-17 . The situation where the two processes spread successively on overlay networks was studied with 
the finding that the outbreak of information diffusion can constrain the epidemic spreading process 18 . An 
analytical approach was developed to provide insights into the symmetric interplay between the two types of 
spreading dynamics on layered networks 19 . A model of competing epidemic spreading over completely overlap- 
ping networks was also proposed and investigated, revealing a coexistence regime in which both types of 
spreading can infect a substantial fraction of the network 20 . 

While the effect of information diffusion (or awareness) on epidemic spreading has attracted much recent 
interest 2128 , many outstanding issues remain. In this paper we address the following three issues. The first 
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concerns the network structures that support the two types of spread- 
ing dynamics, which were assumed to be identical in some existing 
works. However, in reality, the two networks can differ significantly 
in their structures. For example, in a modern society, information is 
often transmitted through electronic communication networks such 
as telephones 29 and the Internet 30 , but disease spreading usually takes 
place on a physical contact network 31 . The whole complex system 
should then be modeled as a double-layer coupled network (overlay 
network or multiplex network) 32 36 , where each layer has a distinct 
internal structure and the interplay between between the two layers 
has diverse characteristics, such as inter-similarity 37 , multiple sup- 
port dependence 38 , and inter degree-degree correlation 39 , etc. The 
second issue is that the effects of one type of spreading dynamics 
on another are typically asymmetric 21 , requiring a modification of 
the symmetric assumption used in a recent work 19 . For example, the 
spread of a disease can result in elevated crisis awareness and thus 
facilitate the spread of the information about the disease 17 , but the 
spread of the information promotes more people to take preventive 
measures and consequently suppresses the epidemic spreading 26 . 
The third issue concerns the timing of the two types of spreading 
dynamics because they usually occur simultaneously on their 
respective layers and affect each other dynamically during the same 
time period 19 . 

Existing works treating the above three issues separately showed 
that each can have some significant effect on the epidemic and 
information spreading dynamics 16,19,40 . However, a unified frame- 
work encompassing the sophisticated consequences of all three issues 
is lacking. The purpose of this paper is to develop an asymmetrically 
interacting spreading-dynamics model to integrate the three issues so 
as to gain deep understanding into the intricate interplay between the 
epidemic and information spreading dynamics. When all three issues 
are taken into account simultaneously, we find that an epidemic 
outbreak on the contact layer can induce an outbreak on the com- 
munication layer, and information spreading can effectively raise the 
epidemic threshold, making the contact layer more resistant to dis- 
ease spreading. When inter-layer correlation exists, the information 
threshold remains unchanged but the epidemic threshold can be 
enhanced, making the contact layer more resilient to epidemic out- 
break. These results are established through analytic theory with 
extensive numerical support. 

Results 

In order to present our main results, we describe our two-layer net- 
work model and the dynamical process on each layer. We first treat 
the case where the double-layer networks are uncorrelated. We then 
incorporate layer-to-layer correlation in our analysis. 

Model of communication-contact double-layer network. Commu- 
nication-contact coupled layered networks are one class of multiplex 
networks 41 . In such a network, an individual (a node) not only 
connects with his/her friends on a physical contact layer (subnet- 
work), but also communicates with them through the (electronic) 
communication layer. The structures of the two layers can in general 
be quite different. For example, an indoor-type of individual has few 
friends in the real world but may have many friends in the cyber 
space, leading to a much higher degree in the communication layer 
than in the physical-contact layer. Generally, the degree-to-degree 
correlation between the two layers cannot be assumed to be strong. 

Our correlated network model of communication-contact layers is 
constructed, as follows. Two subnetworks A and B with the same 
node set are first generated independently, where A and B denote the 
communication and contact layers, respectively. Each layer possesses 
a distinct internal structure, as characterized by measures such as the 
mean degree and degree distribution. Then each node of layer A is 
matched one-to-one with that of layer B according to certain rules. 



In an uncorrelated double-layer network, the degree distribution 
of one layer is completely independent of the distributions of other 
layer. For example, a hub node with a large number of neighbors in 
one layer is not necessarily a hub node in the other layer. In contrast, 
in a correlated double-layer network, the degree distributions of the 
two layers are strongly dependent upon each other. In a perfectly 
correlated double-layer network, hub nodes in one layer must simul- 
taneously be hub nodes in the other layer. Quantitatively, the 
Spearman rank correlation coefficient 39,42 m s , where m s e [ — 1, 1] 
(see definition in Methods), can be used to characterize the degree 
correlation between the two layers. For m s > 0, the greater the cor- 
relation coefficient, the larger degree a pair of counterpart nodes can 
have. For m s < 0, as | m s \ is decreased, a node of larger degree in one 
layer is matched with a node of smaller degree in the other layer. 

Asymmetrically interacting spreading dynamics. The dynamical 
processes of disease and information spreading are typically 
asymmetrically coupled with each other. The dynamics component 
in our model can be described, as follows. In the communication 
layer (layer A), the classic susceptible-infected-recovered (SIR) 
epidemiological model 43 is used to describe the dissemination of 
information about the disease. In the SIR model, each node can be 
in one of the three states: (1) susceptible state (S) in which the 
individual has not received any information about the disease, (2) 
informed state(f), where the individual is aware of disease and is 
capable of transmitting the information to other individuals in the 
same layer, and (3) refractory state (R), in which the individual has 
received the information but is not willing to pass it on to other 
nodes. At each time step, the information can propagate from 
every informed node to all its neighboring nodes. If a neighbor is 
in the susceptible state, it will be informed with probability fi A . At the 
same time, each informed node can enter the recovering phase with 
probability ji A . Once an informed node is recovered, it will remain in 
this state for all subsequent time. A node in layer A will get the 
information about the disease once its counterpart node in layer B 
is infected. As a result, dissemination of the information over layer A 
is facilitated by disease transmission on layer B. 

The spreading dynamics in layer B can be described by the SIRV 
model 26 , in which a fourth sate, the state of vaccination ( V), is intro- 
duced. Mathematically, the SIR component of the spreading 
dynamics is identical to the dynamics on layer A except for different 
infection and recovery rates, denoted by p B and ji B , respectively. If a 
node in layer B is in the susceptible state but its counterpart node in 
layer A is in the infected state, the node in layer B will be vaccinated 
with probability p. Disease transmission in the contact layer can thus 
be suppressed by dissemination of information in the communica- 
tion layer. The two spreading processes and their dynamical inter- 
play are illustrated schematically in Fig. 1. Without loss of generality, 
we set j.i A = fi B = 1. 

Theory of spreading dynamics in uncorrelated double-layer 
networks. Two key quantities in the dynamics of spreading are the 
outbreak threshold and the fraction of infected nodes in the final 
steady state. We develop a theory to predict these quantities for both 
information and epidemic spreading in the double-layer network. In 
particular, we adopt the heterogeneous mean-field theory 44 to 
uncorrelated double-layer networks. 

Let P A {k A ) and P B (k B ) be the degree distributions of layers A and B, 
with mean degree (k A ) and (k B ), respectively. We assume that the 
subnetworks associated with both layers are random with no degree 
correlation. The time evolution of the epidemic spreading is 
described by the variables Sj, A (i), pf A (0> an< ^ 'fc l W> which are the 
densities of the susceptible, informed, and recovered nodes of degree 
k A in layer A at time r, respectively. Similarly, sjt (t)> p% (t), rf B (0> 
and v% (t) respectively denote the susceptible, infected, recovered, 
and vaccinated densities of nodes of degree k B in layer B at time t. 
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(a) Communication-Contact Network (b) t=0 




(c) t=l (d) t=2 

Susceptible Infected Refractory (^^Vaccination 

Figure 1 | Illustration of asymmetrically coupled spreading processes on a simulated communication-contact double-layer network, (a) 

Communication and contact networks, denoted as layer A and layer B, respectively, each of five nodes, (b) At t = 0, node B t in layer B is randomly selected 
as the initial infected node and its counterpart, node Ai in layer A, gains the information that Bi is infected, while all other pairs of nodes, one from layer A 
and another from layer B, are in the susceptible state, (c) At t = 1, within layer A the information is transmitted from A t to A 2 with probability ji A . NodeB 3 
in layer B can be infected by node Bi with probability ft B and, if it is indeed infected, its corresponding node A 3 in layer A gets the information as well. Since, 
by this time, A 2 is already aware of the infection spreading, its counterpart B 2 in layer B is vaccinated, say with probability p. At the same time, node Aj in 
layer A and its counterpart B 1 in layer B enter into the refractory state with probability \i A and respectively, (d) At f = 2, all infected (or informed) 
nodes in both layers can no longer infect others, and start recovering from the infection. In both layers, the spreading dynamics terminate by this time. 



The mean-field rate equations of the information spreading in 
layer A are 



dst (t) 



R A k A ® A (t) + R B ® B (t) Y, k B P B (k B ) 



, (1) 



<(0 A 



dl 



=4M 



p A k A ® A (t) +p B e B (t) Y k B p B (k B ) - P A kA (t),(2) 

The mean-field rate equations of epidemic spreading in layer B are 
given by 



dt 



= -s B kB (t)p B k B ® B (t)-pp A 0 A (t)Y,4 A (mAPA(k A ), (4) 



dp ^=4 B ( t )Pnk B ® B (t)-pl(t), 



Kit) 
dt 



dv? (t) „ . 

=p(J A ® A (t)J24 A (t)kAP A (k A ), 



(5) 

(6) 
(7) 



where ® A (t) (©b(0) is the probability that a neighboring node in 
layer A (layer B) is in the informed (infected) state (See Methods for 
details). 

From Eqs. (l)-(7), the density associated with each distinct state in 
layer A or B is given by 



h,, 



X h (t)=J2 p h(kh)X h kh (t), 



(8) 



where h S {A, B}, X £ {S, J, R, V}, and k h max denotes the largest 
degree of layer h. The final densities of the whole system can be 
obtained by taking the limit t — > °o. 

Due to the complicated interaction between the disease and 
information spreading processes, it is not feasible to derive the exact 
threshold values. We resort to a linear approximation method to get 
the outbreak threshold of information spreading in layer A (see 
Supporting Information for details) as 



where 



P Ac = 



0, for P B >p Bu , 



:{k B )/((kl)-{k B )) 



(9) 



k A 



denote the outbreak threshold of information spreading in layer A 
when it is isolated from layer B, and that of epidemic spreading 
in layer B when the coupling between the two layers is absent, 
respectively. 
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Equation (9) has embedded within it two distinct physical 
mechanisms for information outbreak. The first is the intrinsic 
information spreading process on the isolated layer A without the 
impact of the spreading dynamics from layer B. For fi B > R Bw the 
outbreak of epidemic will make a large number of nodes in layer A 
"infected" with the information, even if on layer A, the information 
itself cannot spread through the population efficiently. In this case, 
the information outbreak has little effect on the epidemic spreading 
in layer B because very few nodes in this layer are vaccinated. We thus 
have ji Bc « p Bu for fl A < f} Au . 

However, for \] A > [i Am epidemic spreading in layer B is restrained 
by information spread, as the informed nodes in layer A tend to make 
their counterpart nodes in layer B vaccinated. Once a node is in the 
vaccination state, it will no longer be infected. In a general sense, 
vaccination can be regarded as a type of "disease," as every node in 
layer B can be in one of the two states: infected or vaccinated. 
Epidemic spreading and vaccination can thus be viewed as a pair 
of competing "diseases" spreading in layer B 20 . As pointed out by 
Karrer and Newman 20 , in the limit of large network size N, the two 
competing diseases can be treated as if they were in fact spreading 
non-concurrently, one after the other. 

Initially, both epidemic and vaccination spreading processes 
exhibit exponential growth (see Supporting Information). We can 
thus obtain the ratio of their growth rates as 



6-- 



(10) 



For 0 > 1, the epidemic disease spreads faster than the vaccination. 
In this case, the vaccination spread is insignificant and can be 
neglected. For 0 < 1, information spreads much faster than the 
disease, in accordance with the situation in a modern society. 
Given that the vaccination and epidemic processes can be treated 
successively and separately, the epidemic outbreak threshold can be 
derived by a bond percolation analysis 20,45 (see details in Supporting 
Information). We obtain 



(*g) 

(i- P s A )({ki)- 



■(k B )Y 



n) 



where S A is the density of the informed population, which can be 
obtained by solving Eqs. (S18) and (S19) in Supporting Information. 
For 0 < 1, we see from Eq. (11) that the threshold for epidemic 
outbreak can be enhanced by the following factors: strong hetero- 
geneity in the communication layer, large information-transmission 
rate, and large vaccination rate. 

Simulation results for uncorrelated networks. We use the standard 
configuration model to generate networks with power-law degree 
distributions 46-48 for the communication subnetwork (layer A). 
The contact subnetwork in layer B is of the Erdos and Renyi (ER) 
type 49 . We use the notation SF-ER to denote the double-layer 
network. The sizes of both layers are set to be N A = N B = 2 X 10 4 
and their average degrees are (k A ) = (k B ) = 8. The degree distribution 
of the communication layer is P A (k A ) = Ck A with the coefficient 



1 



" k A /A and the maximum degree k max ~N 1 ^" /i l \ We 

focus on the case of y A = 3.0 here in the main text (the results for 
other values of the exponent, e.g., y A = 2.7 and 3.5, are similar, which 
are presented in Supporting Information). The degree distribution 

of the contact layer is P B (k B ) = e - <*»> {k B ) k " jk B \. To initiate an 

epidemic spreading process, a node in layer B is randomly infected 
and its counterpart node in layer A is thus in the informed state, too. 
We implement the updating process with parallel dynamics, which is 
widely used in statistical physics 50 (see Sec. S3A in Supporting 
Information for more details). The spreading dynamics terminates 
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Figure 2 | The identification of epidemic threshold on SF-ER networks. 

(a) The susceptibility measure / as a function of the information- 
transmission rate fi A for p = 0.5, ji B = 0.0 (red squares) and fi B = 0. 1 (green 
circles), (b) the threshold fi Ac of information spreading as a function of the 
disease-transmission rate B B for vaccination rate p = 0.5 (red squares) and 
p = 0.9 (green circles), where the red solid lines are analytical predictions 
from Eq. (9). 

when all infected nodes in both layers are recovered, and the final 
densities R A , R B , and V B are then recorded. 

For epidemiological models [e.g., the susceptible-infected-sus- 
ceptible (SIS) and SIR] on networks with a power-law degree distri- 
bution, the finite-size scaling method may not be effective to 
determine the critical point of epidemic dynamics 51,52 , because the 
outbreak threshold depends on network size and it goes to zero in the 
thermodynamic limit 43,53 . Therefore, we employ the susceptibility 
measure 52 % to numerically determine the size-dependent outbreak 
threshold: 



(r 2 )-(r) 2 



(12) 



where N is network size (N = N A = N B ), and r denotes the final 
outbreak ratio such as the final densities R A and R B of the recovered 
nodes in layers A and B, respectively. We use 2 X 10 3 independent 
dynamic realizations on a fixed double-layer network to calculate the 
average value of % for the communication layer for each value of P A . 
As shown in Fig. 2(a), % exhibits a maximum value at ji Ac , which is the 
threshold value of the information spreading process. The simula- 
tions are further implemented using 30 different two-layer network 
realizations to obtain the average value of fi Ac . The identical simu- 
lation setting is used for all subsequent numerical results, unless 
otherwise specified. Figure 2(b) shows the information threshold 
P Ac as a function of the disease-transmission rate ji B . Note that the 
statistical errors are not visible here (same for similar figures in the 
paper), as they are typically vanishingly small. We see that the beha- 
vior of the information threshold can be classified into two classes, as 
predicted by Eq. (9). In particular, for /? B £ ji Bu = ll{k B ) = 0.125, the 
disease transmission on layer B has little impact on the information 
threshold on layer A, as we have P Ac ~P Au = (^a) / {(k 2 ^) — (k A )) 
—0.06. For p B > fS Bu , the outbreak of epidemic on layer B leads to 
P Ac = 0.0. Comparison of the information thresholds for different 
vaccination rates shows that the value of the vaccination probability p 
has essentially no effect on fi Ac . 
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Figure 3 | For SF-ER double-layer networks, epidemic threshold fS Bc as a 
function of the information-transmission rate p A (a) and the vaccination 
rate p (b). In (a), the red solid (p = 0.5) and green dashed (p = 0.9) lines are 
the analytical predictions from Eq. (11), and the blue dot-dashed line 
denotes the case of 0 = 1 from Eq. (10). The inset of (a) shows the 
condition under which a giant residual cluster of layer B exists [from Eq. 
(S20) in Supporting Information] in phase I. In (b), the red solid line (ji A = 
0.05) corresponds to fi Bc = f} Bw and the green dashed line {p A = 0.20) is the 
analytical prediction from Eq. (11). 

Figure 3 shows the effect of the information-transmission rate p A 
and the vaccination rate p on the epidemic threshold P Bc . From 
Fig. 3(a), we see that the value of p Bc is not influenced by \i A for p A 
— Pau ~ 0.06, whereas p Bc increases with P A . For p = 0.5, the 
analytical results from Eq. (11) are consistent with the simulated 
results. However, deviations occur for larger values of p, e.g., p = 
0.9, because the effect of information spreading is over-emphasized 
in cases where the two types of spreading dynamics are treated suc- 
cessively but not simultaneously. The gap between the theoretical 
and simulated thresholds diminishes as the network size is increased, 
validating applicability of the analysis method that, strictly speaking, 
holds only in the thermodynamic limit 20 (see details in Supporting 
Information). Note that a giant residual cluster does not exist in layer 
B forp = 0.9 and p A s 0.49, ruling out epidemic outbreak. The phase 
diagram indicating the possible existence of a giant residual cluster 
[Eq. (S20) in Supporting Information] is shown in the inset of 
Fig. 3(a), where in phase II, there is no such cluster. In Fig. 3(b), a 
large value of p causes p Bc to increase for p A > p Au . We observe that, 
similar to Fig. 3(a), for relatively large values of p, sayp a 0.8, the 
analytical prediction deviates from the numerical results. The effects 
of network size N, exponent y A and SF-SF network structure on the 
information and epidemic thresholds are discussed in detail in 
Supporting Information. 

The final dynamical state of the double-layer spreading system is 
shown in Fig. 4. From Fig. 4(a), we see that the final recovered density 
R A for information increases with p A and p B rapidly for P A £ p Au and 
p B < p Bu . Figure 4(b) reveals that the recovered density R B for disease 
decreases with p A . We see that a large value of p A can prevent the 
outbreak of epidemic for small values of p B , z.&R B ^>0 for p B = 0.2 
and P A a 0.5 (the red solid line). From Fig. 4(c), we see that, with the 
increase in p A , more nodes in layer B are vaccinated. It is interesting 
to note that the vaccinated density V B exhibits a maximum value if p A 
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Figure 4 | For SF-ER networks, the final density in each state versus the 
parameters p A and p B : (a) recovered density R A , (b) recovered density R B , 
(c) the vaccination density V B , and (d) V B versus p B for p A = 0.2, 0.5, 0.9. 

The value of parameter p is 0.5. Different lines are the numerical solutions 
of Eqs. (1) - (8) in the limit f ^> °°, In (a) and (d), we select three different 
values of p A (0.2, 0.5, and 0.9), corresponding to the red solid, green 
dashed, and blue dot-dashed lines, respectively. In (b) and (c), three 
different values of ft B are chosen (0.2, 0.5, and 0.9), corresponding to the 
red solid, green dashed, and blue dot-dashed lines, respectively. 

is not large. Figure 4 shows that the maximum value of V B is about 
0.32, which occurs at p B ~ 0.20, for p A = 0.2. Combining with 
Fig. 3(a), we find that the corresponding point of the maximum value 
P B ~ 0.20 is close to P Bc ~ 0.16 for p = 0.5. This is because the 
transmission of disease has the opposite effects on the vaccinations. 
For p B £ p Bc , the newly infected nodes in layer B will facilitate 
information spreading in layer A, resulting in more vaccinated 
nodes. For p B > P Bc , the epidemic spreading will make a large num- 
ber of nodes infected, reducing the number of nodes that are poten- 
tially to be vaccinated. For relatively large values of p A , information 
tends to spread much faster than the disease for p B £ p Bc , e.g., 0 ~ 
0.21 for p A = 0.5,p = 0.5, p Bc « 0.22, and 0 « 0.12 for P A = 0.9, p = 
0.5, and p Bc ~ 0.23. In this case, the effect of disease transmission on 
information spreading is negligible. The densities of the final dynam- 
ical states for SF-SF networks are also shown in Supporting 
Information, and we observe similar behaviors. 

Spreading dynamics on correlated double-layer networks. In rea- 
listic multiplex networks certain degree of inter-layer correlations is 
expected to exist 35 . For example, in social networks, positive inter- 
layer correlation is more common than negative correlation 54,55 . That 
is, an "important" individual with a large number of links in one 
network layer (e.g., representing one type of social relations) tends to 
have many links in other types of network layers that reflect different 
kinds of social relations. Recent works have shown that inter-layer 
correlation can have a large impact on the percolation properties of 
multiplex networks 37,3 ". Here, we investigate how the correlation 
between the communication and contact layers affects the informa- 
tion and disease spreading dynamics. To be concrete, we focus on the 
effects of positive correlation on the two types of spreading dynamics. 
It is necessary to construct a two-layer correlated network with 
adjustable degree of inter-layer correlation. This can be accomp- 
lished by first generating a two-layer network with the maximal 
positive correlation, where each layer has the same structure as 
uncorrelated networks. Then, Nq pairs of counterpart nodes, in 
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Figure 5 | For two-layer correlated networks with vaccination probability 
p = 0.5, the effect of one type of spreading dynamics on the outbreak 
threshold of the counter-type spreading dynamics, (a) fl Ac versus fi B on 
SF-ER networks with q = 0.0 (red squares) and q = 0.8 (green circles), and 
ER-ER networks with q = 0.0 (blue up triangles) and q = 0.8 (orange down 
triangles). Red solid (SF-ER) and blue dashed (ER-ER) lines are the 
analytical predictions from Eq. (9) and Eq. (S37) (in Supporting 
Information), respectively. The inset shows the inter-layer correlation m s 
as a function of rematching probability q. (b) p Bc versus p A on SF-ER 
networks with q = 0.0 (red squares) and q = 0.8 (green circles), and ER-ER 
networks with q = 0.0 (blue up triangles) and q = 0.8 (orange down 
triangles). Blue solid (q = 0.0) and orange dashed (q = 0.8) lines are the 
analytical predictions for ER-ER networks from Eqs. (S38) - (S41) in 
Supporting Information. 

which q is the rematching probability, are rematched randomly, 
leading to a two-layer network with weaker inter-layer correlation. 
The inter-layer correlation after rematching is given by (see 
Methods) 



m,~\- 



(13) 



which is consistent with the numerical results [e.g., see inset of 
Fig. 5(a) below]. For SF-ER networks with fixed correlation coeffi- 
cient, the mean-field rate equations of the double-layer system 
cannot be written down because the concrete expressions of the 
conditional probabilities P{k B \k A ) and P(k A \k B ) are no longer 
available. 

We investigate how the rematching probability q affects the out- 
break thresholds in both the communication and epidemic layers. As 
shown in Fig. 5, we compare the case of q = 0.8 with that of q = 0.0. 
From Fig. 5(a), we see that q has little impact on the outbreak thresh- 
old p Ac of the communication layer [with further support in Fig. 6(a), 
and analytic explanation using ER-ER correlated layered networks in 
Supporting Information] . We also see that the value of \i Ac for ER-ER 
layered networks with the same mean degree is greater because of the 
homogeneity in the degree distribution of layer A. Figures 5(b) and 
6(b) show that \] Bc decreases with q or, equivalently, ft Bc increases 
with m s . This is because stronger inter-layer correlation can increase 
the probability for nodes with large degrees in layer B to be vacci- 
nated, thus effectively preventing the outbreak of epidemic [see also 
Eqs. (S38)-(S41) in Supporting Information]. Figure 7 shows the 
final densities of different populations, providing the consistent 
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Figure 6 | Effect of varying the rematching probability on outbreak 
thresholds of the two types of spreading dynamics, (a) p Ac versus q on SF- 
ER (red squares) and ER-ER networks (green circles) for p B = 0.05 andp = 
0.5. Red Solid (SF-ER) and green dashed (ER-ER) lines are analytical 
predictions from Eq. (9) and Eq. (S37) in Supporting Information, 
respectively, (b) p Bc versus q on SF-ER (red squares) and ER-ER networks 
(green circles) for p A = 0.2 and p = 0.5. Green solid line is analytical 
prediction for ER-ER networks from Eqs. (S38) - (S41) in Supporting 
Information. 

result that, with the increase (decrease) of q (m s ), the final densities 
R A and R B increase but the density V B decreases. For SF-SF networks, 
we obtain similar results (shown in Supporting Information). 

Discussion 

To summarize, we have proposed an asymmetrically interacting, 
double-layer network model to elucidate the interplay between 
information diffusion and epidemic spreading, where the former 
occurs on one layer (the communication layer) and the latter on 
the counterpart layer. A mean-field based analysis and extensive 
computations reveal an intricate interdependence of two basic 
quantities characterizing the spreading dynamics on both layers: 
the outbreak thresholds and the final fractions of infected nodes. 
In particular, on the communication layer, the outbreak of the 
information about the disease can be triggered not only by its own 
spreading dynamics but also by the the epidemic outbreak on the 
counter-layer. In addition, high disease and information-transmis- 
sion rates can enhance markedly the final density of the informed or 
refractory population. On the layer of physical contact, the epidemic 
threshold can be increased but only if information itself spreads 
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Figure 7 | Effect of rematching probability on the final state, (a) R A versus 
q on SF-ER (red squares) and ER-ER networks (blue up triangles), R B 
versus q on SF-ER (green circles) and ER-ER networks (orange down 
triangles), (b) V B versus q on SF-ER (red squares) and ER-ER networks 
(green circles). Different lines represent the analytic solutions for ER-ER 
networks, calculated by summing the final densities of all degrees from Eqs. 
(S28) - (S34) in Supporting Information. The parameter setting is p A = 
0.2, p B = 0.4 and p = 0.5. 
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through the communication layer at a high rate. The information 
spreading can greatly reduce the final refractory density for the dis- 
ease through vaccination. While a rapid spread of information will 
prompt more nodes in the contact layer to consider immunization, 
the authenticity of the information source must be verified before 
administrating large-scale vaccination. 

We have also studied the effect of inter-layer correlation on the 
spreading dynamics, with the finding that stronger correlation has no 
apparent effect on the information threshold, but it can suppress the 
epidemic spreading through timely immunization of large-degree 
nodes 56 . These results indicate that it is possible to effectively mitigate 
epidemic spreading through information diffusion, e.g., by inform- 
ing the high-centrality hubs about the disease. 

The challenges of studying the intricate interplay between social 
and biological contagions in human populations are generating 
interesting science 57 . In this work, we study asymmetrically interact- 
ing information-disease dynamics theoretically and computation- 
ally, with implications to behavior-disease coupled systems and 
articulation of potential epidemic-control strategies. Our results 
would stimulate further works in the more realistic situation of asym- 
metric interactions. 

During the final writing of this paper, we noted one preprint 
posted online studying the dynamical interplay between awareness 
and epidemic spreading in multiplex networks 58 . In that work, the 
two competing infectious strains are described by two SIS processes. 
The authors find that the epidemic threshold depends on the topo- 
logical structure of the multiplex network and the interrelation with 
the awareness process by using a Markov-chain approach. Our work 
thus provides further understanding and insights into spreading 
dynamics on multi-layer coupled networks. 



contain the variable s^ B ( t) because a node in layer B must be in the susceptible state if 
its counterpart node in layer A is in the susceptible state. 

Spearman rank correlation coefficient. The correlation between the layers can be 
quantified by the Spearman rank correlation coefficient 39,42 defined as 



'n(w 2 -i) 



(16) 



where N is network size and A ; denotes the difference between node is degrees in the 
two layers. When a node in layer A is matched with a random node in layer B, m s is 
approximately zero in the thermodynamic limit. In this case, the double-layer 
network is uncorr elated 39 ' 42 . When every node has the same rank of degree in both 
layers, we have m s ~ 1. In this case, there is a maximally positive inter-layer 
correlation where, for example, the hub node with the highest degree in layer A is 
matched with the largest hub in layer B, and the same holds for the nodes with the 
smallest degree. In the case of maximally negative correlation, the largest hub in one 
layer is matched with a node having the minimal degree in the other layer, so we have 
m s ~ — 1 . 

In a double-layer network with the maximally positive correlation, any pair of 
nodes having the same rank of degree in the respective layers are matched, i.e., A; = 0 
for any pair of nodes A, and B,. We thus have m s = 1, according to Eq. (16). After 
random rematching, a pair of nodes have A, = 0 with probability 1 — q and a random 
difference A', with probability q. Equation (16) can then be rewritten as 



m s — 1 - 



EN a/2 

' N(N 2 -1) 



(17) 



When all nodes are randomly rematched, the layers in the network are completely 
uncorrelated, i.e., m s ~ 0. In this case, we have 



5 w(w 2 -i) 



(18) 



Submitting Eq. (18) into Eq. (17), the inter-layer correlation after rematching is given 
by 

m s -l-q. (19) 



Methods 

Mean-Field theory for the uncorrelated double-layer networks. To derive the 
mean-field rate equations for the density variables, we consider the probabilities that 
node Aj in layer A and node B t in layer B become infected during the small time 
interval [f, t + dt]. On the communication layer, a susceptible node A, of degree k A 
can obtain the information in two ways: from its neighbors in the same layer and from 
its counterpart node in layer B. For the first route, the probability that node A, receives 
information from one of its neighbors is k A fi A ® A (t)dt, where & A {t) is the probability 
that a neighboring node is in the informed state 59 and is given by 



T,K' i (kA-l)P A (kA)pl(t) 



(14) 



where (k A ) — k A P A {k A ). To model the second scenario, we note that, due to the 

asymmetric coupling between the two layers, a node in layer A being susceptible 
requires that its counterpart node in layer B be susceptible, too. A node in the 
communication layer will get the information about the disease once its counterpart 
node in layer B is infected, which occurs with the probability 
1 ^2 l k ^(^l^)^B^B®B( f )^> where P(k B \k A ) denotes the conditional probability that 
a node of degree k A in layer A is linked to a node of degree k B in layer B, and 
k B Rft® B {t)dt is the probability for a counterpart node of degree k B to become infected 
in the time interval [t, t + dt]. If the subnetworks in both layers are not correlated, we 
have P(k B \k A ) = P B (k B ). The mean-field rate equations of the information spreading 
in layer A are Eqs. (l)-(3). 

On layer B, a susceptible node B, of degree k B may become infected or vaccinated in 
the time interval [t, t + dt]. This can occur in two ways. Firstly, it may be infected by a 
neighboring node in the same layer with the probability k B fi B (-) B {t)dU where ® B (i) is 
the probability that a neighbor is in the infected state and is given by 



O s (f) = 



E«,(fe-iWK(t) 
<*■> 



(15) 



where (k B ) — k B P B (k B ). Secondly, if its counterpart node in layer A has already 

received the information from one of its neighbors, it will be vaccinated with prob- 
ability p. The probability for a node in layer B to be vaccinated, taking into account the 
interaction between the two layers, is p P(k A \k B )sf A (t)fi A k A ® A (t)dt, where 
P{k A \k B ) denotes the conditional probability that a node of degree k B in layer B is 
linked to a node of degree k A in layer A, and s£ (t)fi A k A © A (t)dt is the informed 
probability for the counterpart node of degree k A in the susceptible state [P(k A \k B ) — 
Ra(^a) f° r m s — 0]. The mean-field rate equations of epidemic spreading in layer B are 
Eqs. (4) - (7). We note that the second term on the right side of Eq. (4) does not 
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